arXiv:1504.03378v3 [cond-mat.supr-con] 4 Aug 2015 


Random field disorder and charge order driven quantum oscillations in cuprates 
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In the pseudogap regime of the cuprates, charge order breaks a Z 2 symmetry. Therefore, the interaction of 
charge order and quenched disorder due to potential scattering, can, in principle, be treated as a random field 
Ising model. A numerical analysis of the ground state of such a random field Ising model reveals local, glassy 
dynamics in both 2D and 3D. The glassy dynamics are treated as a heat bath which couple to the itinerant 
electrons, leading to an unusual electronic non-Fermi liquid. If the dynamics are strong enough, the electron 
spectral function has no quasiparticle peak and the effective mass diverges at the Fermi surface, precluding 
quantum oscillations. In contrast to charge density, d-density wave order (reflecting staggered circulating cur¬ 
rents) does not directly couple to potential disorder, allowing it to support quantum oscillations. At fourth order 
in Landau theory, there is a term consisting of the square of the d-density wave order parameter, and the square 
of the charge order. This coupling could induce parasitic charge order, which may be weak enough for the Fermi 
liquid behavior to remain uncorrupted. Here, we argue that this distinction must be made clear, as one interprets 
quantum oscillations in cuprates. 

PACS numbers: 74.72.Kf, 7I.10.Hf, 73.22.Gk 


Recent experiments have observed a short-ranged incom¬ 
mensurate charge density wave (ICDW) order in the under¬ 
doped regime of the cuprates uJ-lst], providing a tempting ex¬ 
planation for the underlying order of the pseudogap regime. 
Here, we take a critical view of ICDW as the underlying or¬ 
der in terms of its ability to support quantum oscillations, 
which are generally agreed to reflect a Fermi surface recon¬ 
struction jg Ul in . and therefore a Fermi liquid ground state, at 
least in the sense of continuity ifT^I . To date, there is no gen¬ 
eral agreement as to the precise nature of this reconstruction. 

Because strict ICDW does not have a sharply defined Fermi 
surface 01, there can be no quantum oscillations that are 
truly a periodic function of 1/B (where B is the magnetic 
field). We show that even commensurate charge density wave 
(CDW) order—chosen for simplicity to be of period-2—in 
the presence of disorder may not be able to explain a Fermi 
surface reconstruction and consequently quantum oscillations. 
In short, ubiquitous potential disorder necessarily couples to 
CDW order, leading to a non-Fermi liquid electron spectral 
function without quasiparticles. If this is the case, the princi¬ 
pal order could not be CDW. 

Another possibility for quantum oscillations is the d- 
density wave (DDW) proposed previously llldll . This order, 
illustrated in Fig. [T] in its period-8 version, reflects staggered, 
circulating currents, making it impervious to direct poten¬ 
tial scattering. To the extent that period-8 DDW can in¬ 
duce period-4 CDW, the DDW can be affected by potential 
disorder—^but only at 4th order in Landau theory. Experi- 
the situation is unclear: some neutron scattering 
are consistent with DDW order, but nuclear 
magnetic resonance (NMR) measurements find no circulat¬ 
ing currents iH (see, however, Ref. [l3 for a dissenting opin¬ 
ion). The period-8 DDW has one electron pocket and two 
smaller hole pockets in the reduced Brillouin zone, thus pro¬ 
viding an explanation of the quantum oscillations of the Hall 
coefficients ifl^ . 

For simplicity, we focus on period-2 CDW, which breaks 
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FIG. I. Current pattern for period-8 DDW, reproduced from Ref.[^ 
The wave vector Q = (^, -), where a is the lattice constant I19l 
I20II . In Landau theory, it can couple to CDW of period 2Q. The 
relative magnitudes of the currents are depicted by the thickness of 
the arrows in the legend. Note the antiphase domain wall structure. 


Z 2 symmetry and necessarily couples to disorder. On symme¬ 
try grounds, the effective Hamiltonian is modeled by a random 
field Ising model (RFIM), ^ 

= ( 1 ) 

iU) i 


where J is the coupling between the Ising spins, and {/if} is 
a set of uncorrelated, uniformly distributed (rectangular dis¬ 
tribution) random variables with zero mean and variance 
The notation {ij) denotes nearest neighbors. The model is 
controlled by a single dimensionless parameter, C = 7- 
The disorder in a RFIM drives fluctuations on many length 
scales, and consequently many time scales, producing glassy 
dynamics and a frequency-dependent susceptibility x(a;) El. 
This result is recapitulated here for the 2D case and extended 
to the essential 3D case. Analogous to the thermally driven 
fluctuation of an Ising system at finite temperature, the RFIM 
has disorder driven fluctuations at zero temperature. A dis¬ 
tribution p{L) of domain walls of length L arises from the 
domains in the ground state of the RFIM ll22ll . 
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D xq X 6 ^0 k 

2 1.41(4) 0.28(2) 0.31(2) 0.75(2) 2 

3 0.2(1) 0.33(1) 0.137(6) 0.47(9) 5.6(1) 

TABLE I. Best fit parameters for Pdw in FIG.I^and Eq.[3l 


The appearance of domain walls in RFIM is identified nu¬ 
merically by converting the RFIM to a network flow model 
1 231], solving the “minimum-cut” problem, for which there 
are efficient algorithms 11241] . Briefly, by careful choice of the 
parameters of the flow-network and the addition of two Ac¬ 
tional source and sink nodes, each cut is made to correspond 
to a spin configuration such that the minimal cut corresponds 
to the RFIM ground state configuration. The probability that a 
domain wall of linear dimension L exists in the ground state, 
Pdv/{L), is determined by averaging over many disorder re¬ 
alizations. To help understand the meaning of this quantity, 
notice that, in the ID case, a domain wall is just a spin flip. 
In the disorder-free case, creating a single spin flip costs en¬ 
ergy ^ 2 J, an energy that, by Jordan-Wigner transformation, 
can be thought of as a fermion gap. The spin flip can move 
throughout the system at no energy cost. 

In higher dimensions, the analogy is less precise, but the 
presence of a domain wall results in the collapse of the gap in 
the Ising system. Most importantly, the size of the domains 
is controlled by locations of these domain walls. In particu¬ 
lar, Pdw is the cumulative distribution function of the ordered 
domains or “clusters” of linear dimension L, and therefore 


p{L) = 


dL 


Pdw is found to lie on a universal curve 
metric sigmoid. 


Pdv, ~ f{x) = 


1 


(l + exp[^]) 


( 2 ) 


which is an asym- 


( 3 ) 


where x = logP — (C/Co) ond Co is numerically fit, and sets 
a scale for the strength of the disorder. In 2D, k is set to 2 
as in Ref. in agreement with the analytical result for the 
special case Pdw = 1/2 1^ . In 3D, k is numerically fit. The 
sigmoid’s best fit parameters xo and A control its center and 
width, respectively, while 9 controls the asymmetry. Phys¬ 
ically, Xq determines the onset of the occurrence of domain 
walls, and A how quickly the regime is dominated by the ex¬ 
istence of at least one domain wall. The numerical results are 
summarized in Fig.|2]and Table|I] 

The distribution p{L) is important because it controls the 
(necessarily local) susceptibility |22, 2^ 


Imx(w) ~ J dLp{L)5 — wqc ^ 


( 4 ) 


This phenomenological argument for the susceptibility cap¬ 
tures the essential glassy characteristics resulting from p{L). 
In principle, the attempt frequency loq, the fractal dimen¬ 
sion a, and the length scale c are microscopic parameters. 



logL-(C/Co)2 



FIG. 2. Scaling of Pdw in 2D and 3D. All points are averages over 
2048 disorder realizations, except in 3D for L = 40 and L = 80, 
with 1024 realizations and for L = 160, with 512 realizations. 


which are left undetermined. Notice that the fractal dimen¬ 
sion a < D, where D is the ambient spatial dimension. For 
2D and 3D, the integral simplifies in the small ui limit to 

Imxiuj) ^ — ( 5 ) 

(jJ 

We have put 17 = 1/ (log for clarity and compactness; 
17(w) is strictly increasing for 0 < w < ujq, and vanishes as 
a; —1 0. The exponent 

^ = 1 + 1/ (Aa) > 1 (6) 

depends only on the fractal dimension of the domains a and 
on their distribution of sizes via the parameter A. Moreover, in 
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both 2D and iD, the numerical value of A was found to lead 
to Ip > 2 (see Table|I|i. 

We now focus on the interaction of the itinerant electrons 
with the emergent glassy CDW order, assumed to enter as 
a heat bath of fluctuations of the RFIM. The self energy S 
of the electrons is calculated to leading order in perturbation 
theory (see FIG. O, assuming some coupling 7 of the RFIM 
fluctuations to the electrons, from the form of x in Eq. |5] in 
a reduced graph expansion 10. It is unnecessary to use the 
matrix formalism corresponding to the charge order, because, 
as we shall see, there are no quasiparticles, and hence no pos¬ 
sible Fermi surface reconstruction. In terms of the energy of 
quasiparticles lo, 

/ d ^ 

^ Im G(k — q, w — a;')Im x(q, w') 
^ q 

X IK^') + /(w - w')] (7) 


The Fermi and Bose functions /(w) and b{uj) restrict the w' 
integration to [ 0 , w] in the zero-temperature limit we are con¬ 
sidering, making the integral vanish for w < 0. Because the 
susceptibility is local the self energy is also local, and the sum 
over q reduces to the density of states at the Fermi surface, z/: 


pLO 

ImS(a;) = — 7 ^ 1 / / luixioo') duo' 

^0 


Ip — 1 


where Eq = w > 0. From the Kramers-Kronig 

relations: 


ReE(a;) = -P 


a;'ImE(w') 


duj' 


2 En 


TTplp — 1) Jo ~ 


( 8 ) 


where we have introduced a cutoff A. Because fl is slowly 
varying, we approximate it as a constant with w' = w: 


Re E(w) 


2 So 

TT UJ TTyip — 1) 


where we have taken the largest possible value of the cutoff, 
A —)• Wo, and discarded the subdominant terms in the limit 
w 0. Because '0 > 2, as w —0, both the real and the 
imaginary part of the self energy vanish. 

The spectral function A{k = kp, oj) is plotted in FIG.|4]for 
several values of ip. The emergent behavior is of an unusual 
non-Fermi liquid; for k = kp and in the w —0 limit. 


A{kp,cu)^jt^n^-^. (10) 

4 2^0 

Provided that "0 < 3 or equivalently a > ( 2 A)-i 01 , the 
spectral function vanishes as w —)> 0. The falloff is extremely 
slow, behaving as a fractional power of a logarithm. Further¬ 
more, and despite the slow falloff, the quasiparticle weight 



FIG. 3. Leading order (one-loop) self energy graph. The fermion 
couples to the bath of RFIM fluctuations. 



FIG. 4. Plot of spectral density as a function of frequency relative 
to the chemical potential, with wave vector k = kp- Several values 
of the exponent ip are plotted. The relative scale of each curve is 
arbitrary: the unspecified prefactor Eq is not included. 


always vanishes, and equivalently the effective mass diverges, 
as w ^ 0 : 


Z-i = 1 - Re^ = 1 
duo 


2 Eo Ip _ 

TTW Ip — 1 


( 11 ) 


Cuprates are reasonably modeled as weakly coupled stacks 
of 2D layers 10 . The above work addresses isotropic cou¬ 
pling; we now argue that anisotropy will not materially affect 
the results. Consider the Hamiltonian 


i7stacked = -J\\ Y. ^ ^ 

PAxy PAx * 

where J|| is the in-plane coupling and the interplane cou¬ 
pling. {ij)^ denotes nearest neighbors in the z direction, and 
iU)xy neighbors in the xy plane. The random fields hf 
are as before. 

Unlike 2D, in 3D there is a order-disorder phase transition. 
In the isotropic case, i.e., J = J 2 = J||, the zero temperature 
phase transition occurs at a finite C, = ^, found numerically 
to be Cc = 0.446 ± 0.001, in good agreement with previous 
results 10 . 
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FIG. 5. Phase diagram for the stacked problem. The 3D isotropic 
case corresponds to the dashed diagonal line D = Jy. For each 
value of Jz considered, the numerically identified phase transition 
is a red x. The mean field result, Eq. M is the red line (includ¬ 
ing an approximate scale factor). The ordered phase is above the x 
symbols. 


The anisotropic case (with ^ Jy) is illustrated in Fig. |5] 
Numerically, a particular value of Jz is fixed, and Jy is varied 
to identify the phase boundary in the J^-Jy plane. A sim¬ 
ple mean field theory result is also illustrated: as shown ear¬ 
lier lE^I . in 2D the correlation length 


6d[J/o-] -- exp 



( 13 ) 


with (^0 ~ 0.75. Treating the 3D system as a stack of cou¬ 
pled 2D planes, a mean field theory argument suggests the 
crossover from purely 2D (at weak enough Jz) to 3D occurs 
for 


Jz > J\\/e2D- ( 14 ) 

The qualitative features are readily understood. When 
J|| —0 the system decouples as a ID RFIM, which can¬ 
not order (a scenario irrelevant to the cuprates). On the other 
hand, when Jz ^ 0 the case simplifies to the 2D RFIM, 
which while it also cannot order, has an exponentially large 
crossover scale. It is in the latter regime that the fully 3D and 
stacked 2D results overlap. For weak Jz, the system is dis¬ 
ordered and the total energetic contribution from the Jz cou¬ 
pling can be made small relative to the in-plane Jy terms. An 
interpolation between the 2D and 3D cases is expected in the 
anisotropic case, which should always result in the supression 
of a quasiparticle peak. 

In conclusion, random field disorder is significant even in 
the apparently well ordered materials of high temperature su¬ 
perconductors, but its effect is quite different for the two or¬ 
ders, CDW and DDW. Because the Z 2 symmetry is broken for 


period-2 CDW, it is susceptible to random field disorder, de¬ 
stroying the Fermi surface, as we have found here by treating 
it as a RFIM. The disorder results in the glassy susceptibility 
Imx(w) of Eq. |5] producing a quite unusual non-Fermi liq¬ 
uid. Physically, the glassy dynamics are due to the wide range 
of scales over which domain walls exist in the ground states 
of the 2D and 3D RFIM and are characterized by the parame¬ 
ter Ip, which controls susceptibility and in turn the non-Fermi 
liquid behavior. No Fermi-surface reconstruction can in prin¬ 
ciple occur, precluding quantum oscillations, up to some im¬ 
portant caveats: the coupling parameter Eg must not be too 
small, and the CDW correlation length cannot be too large 
relative to the cyclotron radius (see Ref. m. Truly incom¬ 
mensurate order in the presence of disorder is far too complex 
and was not addressed in the present work, but can only make 
things worse as far as quantum oscillations are concerned. 

In contrast, DDW modulates bond currents—a Hartree- 
Fock calculation of DDW is given by Laughlin ll32n —which 
cannot directly couple to potential disorder, even though the 
order breaks translational symmetry. No non-Fermi liquid be¬ 
havior is expected. Higher periodicity DDW (for example, 
period-8) can induce parasitic charge order that can couple to 
disorder. Being a higher order effect in Landau theory, this 
coupling may be weak. However, the observed weak CDW 
involves such a small motion of the atoms, it is hard to believe 
that it could be the cause of a large magnitude pseudogap. In 
any case, the short range nature of the CDW combined with 
RFIM disorder cannot explain quantum oscillations, at least if 
the resulting electronic state is a non-Fermi liquid. As a third 
option, if we neglect disorder and assume very long-ranged 
CDW, (perhaps infinitely long-ranged), Fermi surface recon¬ 
struction and quantum oscillation have been shown to be pos¬ 
sible [33-3^. The current experiments, however, do not ap¬ 
pear to support such assumptions. 
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